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ABSTRACT 

The high-frequency peaked BL Lac PKS 2155-304 at redshift z = 0.116 is a well-known VHE (>100 
GeV) 7-ray emitter. Since 2002 its VHE flux has been monitored using the H.E.S.S. stereoscopic array 
of imaging atmospheric-Cherenkov telescopes in Namibia. During the July 2006 dark period, the av- 
erage VHE flux was measured to be more than ten times typical values observed from the object. This 
article focuses solely on an extreme 7-ray outburst detected in the early hours of July 28, 2006 (MJD 
53944). The average flux observed during this outburst is I(>200 GeV) = (1.72±0.05 st at±0.34 syst ) x 
10~ 9 cm~ 2 s , corresponding to ~7 times the flux, I(>200 GeV), observed from the Crab Nebula. 
Peak fluxes are measured with one-minute time scale resolution at more than twice this average value. 
Variability is seen up to ~600 s in the Fourier power spectrum, and well-resolved bursts varying on 
time scales of ^200 seconds are observed. There are no strong indications for spectral variability 
within the data. Assuming the emission region has a size comparable to the Schwarzschild radius of 
a ~10 9 M Q black hole, Doppler factors greater than 100 are required to accommodate the observed 
variability time scales. 

Subject headings: Galaxies: active - BL Lacertae objects: Individual: PKS 2155— 304 - Gamma rays: 
observations 
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1. INTRODUCTION 

Flux variability studies provide a strong probe into 
the physical processes of the innermost regions of Ac- 
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tive Galactic Nuclei (AGN). Although the broad-band 
emission from all AGN is highly variable, the most ex- 
treme flux variability, i.e. largest magnitude and short- 
est time scale, is observed from a class of AGN known 
as blazars. As a result blazar variability studies are cru- 
cial to unraveling the mysteries of AGN. Over a dozen 
blazars have been detected so far at VHE energies. In 
the Southern Hemisphere, PKS 2155— 304 is generally 
the brightest blazar at these energies, and is probably 
the bes t-studied at all waveleng ths. The VHE flux ob- 
served (jAharonian et al. 2005af ) from PKS 2155-304 is 
typically of the order ~15% of the Crab Nebula flux 
above 200 GeV. The highest flux previously measured 
in one night is approximately four times this value and 
clear VHE-flux variability has been observed on daily 
time scales. The mos t rapid flux variability m easured for 
this source is 25 min (HE aroman et al. 2005bl ). occurring 
at X-ray energies. The fastest variation published from 
any blazar, at any wavelength, is an event lasting ~800s 
whe re the X-ray flux from Mkn 501 varied by 30% (Xue & 
Cui[20 05) 30 , while at VHE energies doubling timescales 
as fast as ~15 minut es have been observed from Mkn 421 
(|Gaidos et al. 1996[ ). 

The High Energy Stereoscopic System (H.E.S.S.; 
iHinton 2004 ) is used to study VHE 7-ray emission from 
wide variety of astrophysical objects. As part of the nor- 
mal H.E.S.S. observation program the flux from known 
VHE AGN is monitored regularly to search for bright 
flares. During such flares, the unprecedented sensitivity 
of H.E.S.S. (5 standard deviation, a, detection in ^30 s 
for a Crab Nebula flux source at 20° zenith angle) en- 
ables studies of VHE-flux variability on time scales of 
a few tens of seconds. During the July 2006 dark pe- 
riod, the average VHE flux observed by H.E.S.S. from 
PKS 2155—304 was more than ten times its typical value. 
In particular, an extremely bright flare of PKS 2155—304 
was observed in the early hours of July 28, 2006 (MJD 
53944). This article focuses solely on this particular 
flare. The results from other H.E.S.S. observations of 
PKS 2155-304 from 2004 through 2006 will be published 
elsewhere. 

2. RESULTS FROM MJD 53944 

A total of three observation runs (^28 min each) were 
taken on PKS 2155-304 in the early hours 31 of MJD 
53944. These data entirely pass the standard H.E.S.S. 
data-quality selection criteria, yielding an exposure of 
1.32 h live time at a mean z enith angle of 13°. The stan- 
dard H.E.S.S. calibration (I Aharonian et al. 2004] ) and 
analysis tools (Benbow 2005]) are used to extract the re- 
sults shown here. As the observed sig nal is exceptio nally 
strong, the event-selection criteria (Bejibow_2Q_0iJ are 
performed using the loose cuts, instead of the standard 
cuts, yielding an average post-analysis energy threshold 
of 170 GeV. The loose cuts are selected since they have a 
lower energy threshold and higher 7-ray and background 
acceptance. The higher acceptances avoid low-statistics 
issues with estimating the background and significance 
on short time scales, thus simplifying the analysis. The 

30 Xue & Cui I T2005T) also demonstrate that a 60% X-ra y flux 
increase in ~200s observed (Catancsc & Sarnbruna 2000) from 
Mkn 501 is likely an artifact. 

31 The three runs began at 00:35, 01:06 and 01:36 UTC, respec- 
tively. 



on-source data are taken from a circular region of radius 
e cut = 0.2° centered on PKS 2155-304, and the back- 
ground (off-sour ce data) is estimat ed using the Reflected- 
Region method (|Berge et al. 20071) . 

A total of 12480 on-source events and 3296 off-source 
events are measured with an on-off normalization of 
0.215. The observed excess is 11771 events (-2.5Hz), 
corresponding to a significance of 168er calculat ed fo l- 
lowing the method of Equation (17) in Li & Ma (|1983h . 
It should be noted that use of the standard cuts also 
yields a strong excess (6040 events, 159cr) and results 
(i.e. flux, spectrum, variability) consistent with those 
detailed later. 

2.1. Flux Variability 

The average integral flux above 200 GeV ob- 
served from PKS 2155-304 is I(>200 GeV) = 
(1.72±0.05 sta t±0.34 syst ) x 10~ 9 cm" 2 s -1 , equivalent 
to ^7 times t he I(>200 GeV) obsery ed from the Crab 
Nebula (Ioab ; lAharonian et al. 20061 ) . Figure Q] shows 
I(>200 GeV), binned in one-minute intervals, versus 
time. The fluxes in this light curve range from 0.65 I&ab 
to 15.1 Icrab, and their fract ional root mean squ are 
(rms) variability amplitude (jVaughan et al. 20031 ) is 
F V ar = 0.58 ±0.03. This is ^2 times higher than archiva l 
X-ray variability (Zhang et al. 119991 : IZhang et al. 2005h . 
The Fourier power spectrum calculated from Figure [1] 
is shown in Figure [2j The error on the power spectrum 
is the 90% confidence interval estimated from 10 4 
simulated light curves. These curves are generated by 
adding a random constant to each individual flux point, 
where this constant is taken randomly from a Gaussian 
distribution with a dispersion equal to the error of the 
respective point. The average power expected when the 
measurement error do minates is shown as a dashed line 
(see the Appendix in Vaughan et al. 2003f ). There is 
power significantly above the measurement noise level 
up to 1.6 x 10 _3 Hz (600s). The power spectrum also 
shows that most of the power is at low frequencies. 
The grey shaded area shows the 90% confidence level 
obtained by simulating 10 4 light curves with a power-law 
Fourier spectrum P v cx v~ ~ 2 (jTimmer fc Koenig 19951) 
and a random Gaussian error as above. The power 
spectrum derived from the data is thus compatible with 
a light curve generated by a stochastic process with a 
power-law Fourier spectrum of index -2. An index of 
-1 produces too much power at high frequencies and is 
rejected. These power spectra are rem arkab ly similar 
to those derived in X-rays (Zhang et al. 1999) from the 
same source. 

Rapid variability is clearly visible in substructures that 
appear in the light curve, with even shorter rise and de- 
cay time scales than found in the Fourier analysis. In 
order to quantify those time scales, the light curve is 
considered as consisting of a series of bursts, which is 
common for AGN and 7-ray bursts (GRBs) . The "gener- 
alized Gaussian" shape from Norris et al. (|1996l ) is used 
to characterize these bursts, where the burst intensity 
is described by: I(i) = ^4exp[— (\t — i ma x|/o'r,d) K ], where 
£max is the time of the burst's maximum intensity (A); a r 
and <7d are the rise (t < i max ) and decay (t > £ max ) time 
constants, respectively; and k is a measure of the burst's 
sharpness. The rise and decay times, from half to maxi- 
mum amplitude, are T r ^d — [In 2] 1 ^ K a r ^- A peak finding 
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TABLE 1 

The results of the best \ 2 fit of the superposition of five 
bursts and a constant to the data shown in figure [t] the 

CONSTANT TERM IS 0.27 ± 0.03 X 10 -9 cm -2 S _1 (1.1 Icrab)- 



^max 


A 






K 


[min] 


[10~ 9 cm -2 s _1 ] 


N 


N 




41.0 


2.7±0.2 


173±28 


610±129 


1.07±0.20 


58.8 


2.1±0.9 


116±53 


178±146 


1.43±0.83 


71.3 


3.1±0.3 


404±219 


269±158 


1.59±0.42 


79.5 


2.0±0.8 


178±55 


657±268 


2.01±0.87 


88.3 


1.5±0.5 


67±44 


620±75 


2.44±0.41 



tool, using a Markov chain algorithm (jMorhac 2000T ). se- 
lected five significant bursts. A function consisting of 
a superposition of an identical number of bursts plus a 
constant signal was fit 32 to the data. The best fit has a 
X 2 probability of 20% and the fit parameters are shown 
in Table [1] Interestingly, there is a marginal trend for 
k to increase with subsequent bursts, making them less 
sharp, as the flare progresses. The k val ues ar e close to 
the bulk of those found by Norris et al. (1911), but the 
time scales measured here are two orders of magnitude 
larger. 

During both the first two bursts there is clear dou- 
bling of the flux within r r . Such doubling is some- 
times used as a characteristic time scale of flux variabil- 
ity. For compatibility with such estimators, the defini- 
tion of doubling time, T2 = FjAT/AI|, from Zhang et 
al. (|1999f ) was used 33 . Here, AT = Tj - T h AI = lj - F, 
Fj = (Ij + F)/2, with T and I being the time and flux, 
respectively, of any pair of points in the light curve. The 
fastest T2 = 224 ± 60 s is compatible with the fastest sig- 
nificant time scale found by the Fourier transform. Av- 
eraging the five lowest T2 values yields 330 ± 40 s. 

The va riability time scale s of these bursts are among 
(see also lAlbert et al. 20071 ) the fastest ever seen in a 
blazar, at any wavelength, and are almost an order of 
magnitude smaller than previously observed from this 
object. It should be noted that similar time scales 
are found with even smaller binning (e.g. 20 s) of the 
H.E.S.S. light curve, and that many checks of the data 
quality were undertaken to ensure that the flux varia- 
tions cannot be the result of background fluctuations, 
atmospheric events, etc. In addition, all the results have 
been verified using an independent calibration method 
and alternative analysis techniques. 

2.2. Spectral Analysis 

Figure [3] shows the time- averaged photon spectrum for 
these data. The data are well fit, x 2 = 17.1 for 13 degrees 
of freedom (d.o.f.), by a broken power-law function: 

E < E ® '■ m = Io (itIv, 
h 



-ri 



p ^ b 1 • dN _ t 1 e b ^^-ri) 1 e \ 

B • \TTev) ITTcvJ 

where I Q = (2.06±0.16±0.41) x 10" 
Eb — 430i22i80 GeV, r x =2.71±0.06±0.10, and 
T2 = 3.53±0.05±0.10. For each parameter, the two 
uncertainties are the statistical and systematic val- 
ues, respectively. Fits to the data of either a simple 

32 The Markov chain burst positions were used to initialize t max 
for each burst. All parameters are left free in the fit. 

33 Only values of T2 with less than 30% uncertainty are consid- 
ered. 



power law (r = 3.19±0.02±0.10, \ 2 = 138, 15 d.o.f) 
or a power law with an exponential cut-off (x 2 = 45, 
14 d.o.f.) are not acceptable. The time-averaged 
spectrum (r = 3.32) of P KS 2155-304 measured in 2003 
(|Aharonian et al. 2005a| ). multiplied by the ratio (48.7) 
of I(>200 GeV) from the respective data sets, is also 
shown in Figure [3] Despite a factor of ~50 change in 
flux there is qualitatively little difference between the 
two spectra. Indeed, fitting a broken power law to the 
current data set, keeping Ti and T2 fixed to the values 
measured in 2003, yields a value for E-q consistent with 
that measured in 2003. The small difference is surprising 
since a change of the spectral shape with varying flux 
levels, typically hardening with increased flux, has often 
been observed from blazars at X-ray energies (see, e.g., 
iGiommi et al. 19901). as we ll as in the VHE domain (see, 
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e.g. 

The high flux observed from PKS 2155-304 allows the 
determination of accurate photon spectra on time scales 
of the order of minutes. Therefore, a simple search for 
temporal changes of the VHE spectral shape within these 
data was performed. Spectra were determined for con- 
secutive data slices of 28 minutes (1 run), 10 minutes, 
and 5 minutes. Fitting the time-average spectral shape, 
allowing only the normalization (J ) to vary, to these 
short-time-scale spectra yields reasonable \ 2 probabili- 
ties. Thus, there are no strong indications of fast spec- 
tral variability. However, weak variations (Ar < 0.2) are 
not ruled out. A more sophisticated study of any fast 
spectral variations within these data is beyond the scope 
of this letter and will be published elsewhere. 

3. DISCUSSION 

It is very likely that the electromagnetic emission in 
blazars is generated in jets that are beamed and Doppler- 
boosted toward the observer. Super luminal expansions 
observed with VLBI (jPiner &: Edwards 2004f ) provide ev- 
idence for moderate Doppler boosting in PKS 2155—304. 
Causality implies that 7-ray variability on a time scale 
tvarj with a Doppler factor 34 (6), is related to the radius 
(R) of the emission zone by R < ct vea 6/ (1 + z). Conser- 
vatively using the best-determined rise time (i.e. r r with 
the smallest error) from Table [T] for i var = 173 ± 28 s 
(note that this is similar to the fastest T2) limits the size 
of the emission region to RS -1 < 4.65 x 10 12 cm < 0.31 
AU. 

The jets of blazars are believed to be powered by accre- 
tion onto a supermassive black hole (SMBH). Thus ac- 
cretion/ejection properties are usually presumed to scale 
with the Schwarzschild radius Rs of the SMBH, where 
Rs = 2GM/c 2 , which is the smallest, most- n atural size 
of the system (see, e.g..[Blandfor d et al. 19 82). Express- 
ing the size R of the 7-ray emitting region in terms of 
Rs, the variability time scale limits its mass by M < 
(c 3 t var <5/2G(l + z))R s /R ~ 1.6 x 10 7 M Q 5R S /R. The re- 
p orted host galaxy lum inosity Mr = —24.4 (Table 3 
in iKotilainen et al. 19981) would imply a SMBH mass of 
order l-2x!0 9 M g (|Bettoni et al. 2003h . and therefore, 
S > 60 — 120i?/i?s- Emission regions of only a few Rs 

34 With <5 defined in the standard way as [T(l— /3 cos 9)]~ 1 , where 
r is the bulk Lorentz factor of the plasma in the jet, /3 = v/c, and 
9 is the angle to the line of sight. 

35 See Wurtz et al. fl996h and Mr > -23.1 (for h = 0.5) showing 
the need for confirmation of this value. 
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would require values of 5 much greater than those typi- 
cally derived for blazars (<5~10) and come close to those 
used for GRBs, which would be a challenge to under- 
stand. For example, the sub-parsec VHE 7-ray emitting 
plasma would have to decelerate with a high efficiency to 
accommodate re latively small Lorentz f actors observed 
at parsec scales (|Piner fc Edwards 2004h . It is however 
possible that the SMBH mass is over-estimated, reducing 
the S constraint by the same factor, or that the variability 
has an origi n (e.g., a geometric e ffect from jet bending as 
discussed in IWagner et al. 1 .9931 unrelated to the black 
hole. Detailed modeling of the spectral energy distribu- 
tion of PKS 2155-304, during the multiple VHE flares 
observed by H.E.S.S. in the July 2006 dark period, in- 
cluding simultaneous multi-frequency data, will appear 
elsewhere. 

The VHE variability observed in this particular flar- 
ing episode is the fastest ever observed from a blazar. 
While the variability is a factor of five times faster tha n 
previously measur ed from Mkn 421 (jGaidos et al. 1 996). 
in terms of the light-crossing time of the Schwarzschild 
radius, Rs/c, the variability of PKS 2155— 304 is another 
factor of « 6 —12 more constrai ning assuming a 10 8 22 M® 
for Mkn 421 (|Woo et al. 20051 ). It should also be noted 
that the choice of a ^3 minute variability time scale here 
is conservative and that the light curve is strongly over- 
sampled, allowing for the first time in the VHE regime 



a detailed statistical analysis of a flare, which shows re- 
markable similarity to other longer duration events at 
X-ray energies. From such rapid variability one must 
conclude that either very large Doppler factors can be 
present in AGN jets, or that the observed variability is 
not connected to the central black hole, clearly showing 
the power of Cherenkov-telescope arrays in probing the 
internal mechanisms in BL Lacs. 
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Fig. 1. — The integral flux above 200 GeV observed from PKS 2155— 304 on MJD 53944 versus time. The data are binned in 1-minute 
intervals. The horizontal line represents I(>200 GeV) observed (Aharonian ct al. 2006) from the Crab Nebula. The curve is the fit to these 
data of the superposition of five bursts (see text) and a constant flux. 
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Fig. 2. — The Fourier power spectrum of the light curve and associated measurement error. The grey shaded area corresponds to the 
90% confidence interval for a light curve with a power-law Fourier spectrum P v <x v~' 2 . The horizontal line is the average noise level (see 
text) . 
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Fig. 3. — The time-averaged spectrum observed from PKS 2155—304 on MJD 53944. The dashed line is the best y 2 fit of a broken power 
law to the data. The solid line represents the fit to th e time-averaged spectrum of PKS 2155—304 from 2003 ( Aharonian et al. 2005a) 
scaled by 48.7. Neither spectrum is corrected (see, e.g., Aharonian et al. 2005b) for the absorption of VHE 7-rays on the Extragalactic 
Background Light. 



